plot[ {sqrt[x (x-1) (x+1)]==y, -sqrt[x (x-1) (x+1)]==y}, {x,-2,2},{y,-2,2} ]